Phase transitions and soft elasticity of smectic elastomers 
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Smectic-C elastomers can be prepared by crosslinking, e.g., liquid crystal polymers, in the smectic- yl 
phase followed by a cooling through the smectic-A to smectic-C phase transition. This transition 
from Daoh to C2h symmetry spontaneously breaks rotational symmetry in the smectic plane as does 
the transition from a smectic-A to a biaxial smectic phase with D2h symmetry. We study these 
transitions and the emergent elasticity of the smectic-C and biaxial phases in three related models 
and show that these phases exhibit soft elasticity analogous to that of nematic elastomers. 

PACS numbers; 83.80.Va, 61.30.-v, 42.70.Df 



Liquid crystalline elastomers j]| are remarkable materi- 
als that combine the orientational and positional order 
of liquid crystals (2(] with the elastic properties of rubber. 
The traditional liquid crystalline nematic, cholesteric, 
smectic- A (SmAV, smectic-C (SmC), and smectic-C* 
phases all exist [l| in elastomeric forms. In this letter 
we investigate the properties of SmC and biaxial phases 
in elastomers formed via spontaneous symmetry break- 
ing from Sm^ or uniaxial phases and, within mean-field 
theory, the phase transitions to them. We introduce and 
analyze phenomenological models for these transitions in- 
volving strains only and a model involving strains and the 
Frank director specifying the direction of local molecular 
order. Our primary result is that monodomain samples 
of the emerging biaxial and SmC phases (see Fig. 0, 
both of which break the continuous rotational symmetry 
in the smectic layers, exhibit soft elasticity characterized 
by the vanishing of certain elastic moduli and the associ- 
ated absence of restoring forces to strains along specific 
symmetry directions. As in monodomain nematic elas- 
tomers ^ IBl this soft elasticity is a consequence of the 
Goldstone theorem that requires any phase with a spon- 
taneously broken continuous symmetry to have modes 
whose energy vanishes with wavenumber. 

To keep our discussion as simple as possible, we 
will consider only smectic elastomers crosslinked in the 
Sm^ phase so that the smectic layers are locked to the 
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FIG. 1: Sample distortion and rotations of the Frank director 
n, the uniaxial anisotropy axis e, and the layer normal N in 
a transition from a (a) SmA to (b) a sheared SmC elastomer. 
In this geometry, the the smectic layers do not rotate. The 
figure corresponds to a shear strain 7 « ul^ < (see text). 



crosslinked matrix p. Macroscopically, these SmA elas- 
tomers are simply uniaxial rubbers (or solids) with Dooh 
symmetry, and we will treat them as such, though we will 
distinguish between the normal N to smectic layers and 
the direction of uniaxial anisotropy e. 

We employ the usual Langrangian formalism 7] in 
which mass points in the undistorted medium, which 
we take as the reference space, are labelled by vectors 
X. Mass points of the distorted medium are at posi- 
tions R(x) = X + u(x) in physical space, which we call 
the target space. Distortions of the reference medium 
are described by the Cauchy deformation tensor A with 
components Ay = dRi/dxj = djRi, i,j — x,y,z. La- 
grangian elastic energies are expressed in terms of the 
nonlinear strain tensor u with components Uij{x.) = 

\{Kjf.Kkj ~ 5ij) = ^{diUj+djUi + diUkdjUk), which trans- 
forms as a tensor in the reference space and which is in- 
variant with respect rotations in the target space. The 
elastic free energy density of a uniaxial elastomer to har- 
monic order in strains can be expressed as 
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where we use the Einstein convention on repeated indices 
and where z is along the uniaxial axis and indices at the 
beginning of the alphabet a,b, .. run over x and y only. 
The tensor Uab — Uab — ^SatUcc is the two-dimensional 
symmetric, traceless strain tensor. The elastic constant 
Ci describes dilation or compression along z. C4 and 
C5 respectively describe shears in the plane perpendicu- 
lar to the anisotropy axis and in the planes containing it. 
Relative volume change SV/V — det A — 1 can be approx- 
imated by Uii in our Landau expansion in powers of , 
and the incompressible limit relevant to most elastomers 
corresponds to C3 — > 00. 

The harmonic free energy /uni describes any uniaxial 
solid or elastomer, including Sm^ elastomers. However, 
to provide a complete description of SmA and SmC elas- 
tomers, we need to add terms that describe the smectic 
layers and the Frank director and their interactions with 
each other and with the elastic medium. Before adding 
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these terms, we will first investigate phase transitions in 
uniaxial elastomers in which either of the elastic con- 
stants C4 or C5 go to zero. The vanishing of C4 and 
C5 lead, respectively, to phases with (orthorhombic) 
and C2h (triclinic) symmetry. We will then consider the 
smectic layers and director and show that the transition 
to the SmC phase, which has C2h symmetry, is in fact 
identical to the one in which C5 goes to zero. 

Strain-only theory for C4 — * 0: If C4 or C5 become 
negative, as they will in response to an instability toward 
biaxial or SmC ordering of constituent mesogens, terms 
higher order in strains must be added to Eq. to ensure 
mechanical stability. If C4 becomes negative, order of the 
shear strain Uab sets in and higher order terms featuring 
Uab have to be added which leads to the model elastic 
energy 

f^ni = /uni + Ai u,,ul^ + A2 u,,ul^, + B {ul^,f, (2) 

where we have dropped qualitatively inconsequential 
higher order terms. It will be useful in the analysis that 
follows to regroup the terms in /^^| by completing the 
squares in ^Ciuj.^ + AiUzzU^i,, etc., and to reexpress it 
as a sum of two terms 

/u!.i'^ = 5 Civl + C2 vzzvu + i C3 + C, ul^ (3) 
fl^^^ =CA + Br{uI,)\ (4) 

where v^z = u^z - aul^, vu = uu - piil^, and where 
/3, which vanishes as C3 00, Br and a, which do 
not, are combinations of the coefficients in fii|. For a 

given iiab, /uni minimized when Uzz = 0iu1,b s-nd uu = 
/3u^j. The equilibrium value of Wab, w°f, = S{caCb - ^Sab), 
where c is a unit vector in the xy plane, is determined by 
minimizing f^^^ '. The result is S* = for C4 > 0, and 
S = ±y/ —C4/ Br for C4 < 0. The equilibrium strain 
Vi^ = \ ( A°"^ A" — 5) is diagonal with components u° = 
'\s+\ij3-l,)S^~uly = -Is + \{I3 - a)S^ , and u^, = 
^aS^ . The new state is biaxial with D2h symmetry. 

To determine the elastic properties of the new state, 
we expand /^^j in powers oi 5u — u — iP . Since the equi- 
librium values oivzz, vn and Uaz are zero, the expansion 
of f.^if ■* is trivial. The structure of f'if is identical to 

•' uni •' uni 

that of an xy model, and it has no restoring force in the 
ordered phase for 5u^y: Sfj;^];P = BrS'^{5u^x - SuyyY . 
Thus it is clear that <5/^^j does not depend on 5uxy to 
harmonic order, i.e., the system is soft with respect to 
shears in the xy plane of the original reference material. 

The strain 5u describes distortions relative to the 
new biaxial reference state measured in the coordinates 
of the original uniaxial state. It is customary, how- 
ever, to express the elastic energy in terms of a strain 
u' — (A°"^)~^fo (A°)~^ measured in the coordinates 
x^ = Xi + = A'^jXj of the new state. A° is not uniquely 
determined by vP; rotations in the target space change 



A but do not change u". In the present case, it is nat- 
ural to choose a" to be diagonal, i.e., not to rotate the 
strain strain after the transition. The elastic energy of 
the biaxial state to harmonic order in the new strains is 
then 

fOih ~ 5 ^zzzz {u'zz)"^ + C^xzxz (u'j.z)^ + ^yzyz {u'y^Y' 
+ Czzxx Uzz'^^xx ~^ ^zzyy U^^Uyy -\- ^ Cxxxx {Uxx) 

~^ 2 ^yyyy ^'^yy^ ~^ ^'^^yy^xx'^yy^ (^) 

where the elastic constants depend on the original elastic 
constants featured in Eq. ||2Jl and the order parameter 
S, and for S" — ^ this energy reduces to the uniaxial 
energy If we take the incompressible limit, C3 — > 00, 
the specifics of the new elastic constants are affected but 
the form of the elastic energy © remains the same. 

Because there was no Su^y term in the expansion of 
/unii there is no term proportional to u'^y as there would 
be in conventional orthorhombic systems. Thus, there is 
no restoring force to xy-stresses, i.e., to opposing forces 
along ±x applied to opposite surfaces perpendicular to y 
or opposing forces along ±y applied to opposite surfaces 
perpendicular to x. In addition, as is the case for stresses 
perpendicular to the anisotropy axis in a nematic elas- 
tomer 0, J it requires no stress to stretch the sample 
along y direction, up to a critical strain value. The same 
soft mode and many more were predicted by Warner and 
Kutter 8] for biaxial nematics forming spontaneously 
from an isotropic elastomer. 

Strain- only theory for C5 — *■ 0: When C5 is driven neg- 
ative the uniaxial state becomes unstable to shear in the 
planes containing the anisotropy axis, and the uniaxial 
energy ^ must be augmented with higher order terms 
involving Uaz to stabilize the system, 

/uni = /uni + ^zzul^ -\- D2 Uuul^ + D3 UabUazUbz 

+ E{ulf, (6) 

where we omitted all unimportant symmetry-compatible 
higher order terms. To study the ordered phase of this 
free energy when C5 < 0, we proceed in much the same 
way as we did for the biaxial state of fj^^. We complete 
squares to write /^^; as the sum of two terms: 

/(^'i^' = i Ci wl + C2 wuwzz + i C3 wl + C4 wl, (7) 

/i^nf^ =C5uL+i?fl("L)', (8) 

where Wzz = Uzz ~ crul^, wu = uu - r-u^^, and Wab = 
iiab - t^iuazUbz ~ ^^abulj whcrc cr, T as C3 
00), uj, and Er are functions of the parameters in f^^^. 
The equilibrium value of Uaz is determined by minimizing 
/uni ' which has xy symmetry. For C5 > 0, Uaz = 0; for 
C5 < and order along x, = S = ±yJ~Cr^/{2ER), 
and M° = 0. The other components of yP are = crS"^ , 
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"^xx = ^{t + uj-<t)S'^, and Uyy = ^{t - uj - a)S'^ . Unlike 
the biaxial case, vP is not diagonal; it has nonvanishing 
xz and zx components that lead to C2h rather than I?2/i 
symmetry. Expanding /^^| in powers of Su = u — vP, 

we find that Sf^^^P ~ AEjiS'^{5uxzY is independent of 
5uyz^ and we might naively expect the system to exhibit 
softness with respect to Uyz- This, however, is not the 
case because /unf depends on 5uyz via the C4 term: 
2C4{duxy — uSduzyY- Thus, the softness of the ordered 
phase with C2h symmetry is more subtle than that of the 
biaxial phase with D2h symmetry. 

To determine the energy of strains relative to the new 
ground state, we need to choose how we define our new 
coordinate system relative to it. It is easiest to visu- 
alize the new state as emerging from a simple shear as 
shown in Fig. ^in which A^^ = dR^/dz = 7 is nonzero 
but hP^x = dR^/dx — 0. In this case, the only nonzero 
components of A" are the diagonal components and A"^ 
which can be expressed in terms of yP: ~ ■\/l + 2u^, 

Kv = \/l + 2<' = 7 = 2<z/Kx, and A^, = 
+ 21*2^ — 7^. The elastic energy can now be written 
in terms of ^' = (A°^)"M^ (A°)"^ as 

fct = 5 C [cos 9 u'^y + sin e J ' + i Czzzz {<z? 

+ Cxzxz {"^xz) "I" Czzxx U^z'^xx + Czzyy U^^Uyy 

+ 2 Cxxxx i^xx) ^" 2 ^yyyy (^yy^ Cxxyy '^xx'^yy 

+ CxxXZ Uxx'^XZ "I" Cyyxz UyyUxz + CzzXZ '^zz'^XZ^ i^) 

where the angle 9 and the elastic constants C, Czzzz and 
so on depend on the original elastic constants in Eq. © 
and S so that one retrieves the uniaxial energy for 
5^0. To lowest order in S, tan 6' = ljS. The elas- 
tic energy © has only 12 (including 9) rather than the 
13 elastic constants of conventional triclinic solids 0. 
There are only two rather than three independent elas- 
tic constants in the subspace spanned by u'^y and u'y^. 
With the introduction of w = {u'^y,u'^y), elastic energies 
in this subspace can be expressed as ■ m ■ v where m 
is a 2 X 2 matrix. In general m has two independent 
eigenvalues mi and m2 with respective associated orth- 
normal eigenvectors ei and 62 with respect to which m 
is diagonal: v ■ m ■ v — mi{ei ■ vY + "^2(62 • ^Y- The 
first terms in Eq. is of the form \mi{ei ■ vY with 
mi — C and ei = (cos 0, sin 0), and we conclude that 
TO2 — 0. Thus distortions along e?2 — {— sa\9,cos9), i.e., 
distortions for which v \\ 6*2 cost no energy (see Fig. 
Stated differently, there are no restoring forces to stress 
— sm9(Txy + cos 9 a zy — e2i<Jiy, i.e., to stress in the xz 
plane directed along 62 or stress in the plane perpendic- 
ular to 6*2 directed along y. 

Theory with director and smectic layers: The following 
theory generalizes the achiral limit of a continuum the- 
ory for SmC* elastomers by Terentjev and Warner [loj 



in a formalism that ensures invariance with respect to 
arbitrary rather than infinitesimal rotations of both the 
director and mass points. 

In traditional uncrosslinked liquid crystals, there is no 
reference space, and all physical fields like the smectic 
layer-displacement field U, the layer normal N, and the 
Frank director n are defined at real or target space points 
R and they transform as scalars, vectors, and tensors 
under rotations in the target space. In the Lagrangian 
theory of elasticity, fields are defined at reference space 
points X, and they transform into themselves under the 
symmetry operations of that space. To develop a theory 
of liquid-crystalline elastomers, it is necessary to com- 
bine target-space liquid crystalline fields and reference 
space elastic variables to produce scalars that are invari- 
ant under arbitrary rotations in the target space and un- 
der symmetry operations of the reference space. This 
requires that we be able to represent vectors and ten- 
sors in either space 1^. The matrix polar decomposition 
theorem ^l| applied to the Cauchy tensor A provides a 
route to this representation. Like any non-singular ma- 
trix, A can be decomposed into the product of an orthog- 
onal matrix O times a symmetric matrix: A = OM}^'^, 

1 /2 

where M_ ' is the symmetric square root of the sym- 
metric matrix M = ( A"^ A) = (1 -I- 2u) , which depends 

_ _ 1/2 

only on the symmetric strain u, and O = AA£ ■ The 

orthogonal matrix O converts (or rotates) any reference 
space vector a to a target space vector a via a = O ■ a 
and a target-space vector to a reference space vector via 
a = ■ a. Equipped with O, we can rotate the target 
space director n to a reference space vector n = (c,n^), 
where riz — a/I — c^, and rotate the reference space 
anisotropy vector e = (0,0,1) to a target-space vector 
e = O • e. From these, we can form invariant couplings 
like (fi • e)^ = (n • e)^ = 1 - c^, CaUabCb, and CaUazfiz as 
building blocks the elastic energy. 

To treat smectic layers, we need to discuss in more 
detail the smectic displacement field U and the layer 
normal N. The smectic mass-density-wave amplitude 
for a system with layer spacing d has a phase (^(R) = 
<lo[Rz — C^(R-)] where = 27r/d. Since there is a one- 
to-one mapping from the reference space points x to 
the target-space points R(x), we can express as a 
function of x as 0(x) = qQ[z + Uz{x) — ?7(R(x))]. We 
are only considering systems crosslinked in the smectic 
phase. In these systems, the smectic mass-density wave 
cannot translate freely relative to the reference material, 
and there is a term in the free-energy density \ A{uz — UY 
that locks the smectic field U to the displacement field 
Uz . In what follows, we will take this lock- in as given 
and set U — Uz- The layer normal in the target space 
is N, = V,0/|V0|, where = d(j)/dRi = Sj^A"/. 
Thus, when U is locked to Uz, Ni = [(M"^)^z]"^/^Aj/, 
Ni = [{M.~^)zz]~^^'^ M^P^^ , and to harmonic order in 
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and , N ■ n = N • n = 1 - - u^aCa H • 

We can now develop a full phenomenological free en- 
ergy for the Sm^-to-SmC transition in an elastomer. 
In the equilibrium Sm^ phase, the director is parallel to 
both the layer normal N and the anisotropy axis e, which 
are parallel to each other, and there are energy costs pro- 
portional to (N-n)^ and (e-n)^ associated with deviations 
from this equilibrium. These combine to yield a term 
in the free-energy density proportional to and higher 
order terms involving the strain and strain-director cou- 
pling. With the addition of higher order terms in (N ■ n)^ 
and (e • fi)^ to stabilize the SmC phase, the free energy 
up to inconsequential higher order terms becomes 

/ = /uni + /tilt + /coupl, (10) 

where /uni is the uniaxial energy and 

hit = hrcl + \g{il)\ (11) 

/coupl = Ai C^Uz2 + A2 C^Ubb -I- A3 CaUabh 

+ \4CaUazflz- (12) 

When Ai = A2 = A3 = 0, this model is equivalent to that 
studied in Ref. jlOf when polarization is ignored. 

We can now analyze the transition to the SmC phase 
in exactly the same way as we did in the strain only 
model. We complete the squares involving the strains 
and the director-strain couplings. The result is an ex- 
pression for the free energy that is the sum of a term 
/^^^ quadratic in shifted strains of the form Uzz + pc^, 
etc. and a term /^^^ — ^rucl + ^ gn (c^)^ depending 
only on c^, where rjf — r — Xl/{2Cii). With c chosen to 
align along x, minimization yields c'^ = S = —rnj gn, 
= —XiS/ and diagonal components of ■if' pro- 
portional to S"^ . As anticipated from the strain-only 
model for C5 0, the SmC phase is a sheared biax- 
ial phase with C2h symmetry. If Ca is integrated out of 
/, the result is identical to /uni with C5 renormalized to 
Cs.iJ = C's — A4/(2r), which vanishes at = 0. 

To address the elastic properties of the SmC phase 
within our last model, we expand the complete elastic 
energy IjlUI) in terms oi 5u = u — iP and 5c = c — c°. 
Since /^^^ has xy symmetry, its expansion has a [Scx)'^ 
but no {6cyY term. The expansion of /^^^ has terms 
quadratic in 5u and couplings between 5u and 5ca- Since 
5cx is massive, it can be integrated out. Converting the 
resulting expression to the strain variable u! of the SmC 
state, we obtain an elastic energy that is identical in form 
to Eq. ((SJ but with an additional term proportional to 
[5cy + p{u'xy + C,u'y^)Y , where ^ and Q are dimensionless 
coefficients that depend on the original elastic constants 
and S. This expression shows clearly that 5cy can relax 
locally to — ^ {u'^y + C,u'y^ to eliminate the dependence 
of the elastic energy on u'^y + C,u'y^ and produce an energy 
with softness identical to that of Eq. @ . 

The above analysis provides us with the strain and 
the c-director c*^ in the coordinates of the reference space 



without rotation. To visualize what happens to e, N, 
and n in the SmC phase, it is useful to consider the 
simple shear of Fig. ^ iii which tan = AJ^^ = 7, and 
A°^ = 0. To lowest order in 5, 7 = 2ul^ = -X^^S/C^, 
and Oij = 5ij -\- ^eizjj. Then n « {nx,Q,l) where 
rix = S + ^7, and = Oiz = (^7,0,1). Thus, at 
the transition e and n rotate relative to N by different 
amounts. The layer normal vector N, however, remains 
fixed in the chosen geometry: Ni ~ A~/ only has a non- 
vanishing z component and is equal to (0,0,1). Thus, 
the smectic layers remain parallel to the a:-axis as one 
would intuitively expect in the chosen geometry. 

Concluding remarks: We have presented models for 
transitions from uniaxial SmA elastomers to biaxial and 
SmC elastomers, and we have calculated the nature of 
the soft elasticity, required by symmetry, of monodomain 
samples of these phases. We hope that our work encour- 
ages experiments to probe this soft elasticity. 
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